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A three-level A system in Tm"^"'' doped YAG crystal is experimentally investigated in the prospect 
of quantum information processing. Zeeman effect is used to lift the nuclear spin degeneracy of 
this ion. In a previous paper [de Seze et al. Phys. Rev. B, 73, 85112 (2006)] we measured 
the gyromagnetic tensor components and concluded that adequate magnetic field orientation could 
optimize the optical connection of both ground state sublevels to each one of the excited state 
sublevels, thus generating A systems. Here we report on the direct measurement of the transition 
probability ratio along the two legs of the Lambda. Measurement techniques combine frequency 
selective optical pumping with optical nutation or photon echo processes. 

PACS numbers: 42.50.Gy, 42.50.Md, 71.70.Ej 



I. INTRODUCTION 

There is an ongoing interest for macroscopic quantum 
processes, especially in the perspective of quantum mem- 
ories and several groups are currently working on per- 
forming storage and retrieval of a quantum state of light 
into an atomic system. Most quantum memory pro- 
cesses rely on the storage of a quantum signal into an 
atomic spin coherence, which is free from spontaneous 
emission. The optical excitation carrying the informa- 
tion resonantly excites an atomic transition. A control 
pulse then converts the optical coherence into a hyperfine 
long-lived coherence. Another control pulse can eventu- 
ally change the spin coherence back into an optical co- 
herence and allow the retrieval of the stored signal. This 
quantum memory scheme can be achieved in a three-level 
A system, where two spin states are connected to a com- 
mon upper level by optical transitions (cf Fig.[T|). 

The mapping of a quantum state of light over a macro- 
scopic atomic ensemble has been demonstrated in various 
materials such as atomic vaporsi or cold atom clouds^. 
The full quantum memory scheme has recently been 
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FIG. 1: Three-level A system coupled by two lasers, fii and 
0.2 are the Rabi frequencies that characterize the atom-laser 
interaction. 



achieved in cold atomic clouds of rubidium''. However, 
rare-earth ion-doped crystals also appear as promising 
candidates for quantum storage applications. They offer 
properties similar to atomic vapors with the advantage 
of no atomic diffusion. At low temperature (< 4 K), the 
optical coherence lifetime can reach several ms, and a hy- 
perfine coherence lifetime has been extended to tens of 
seconds in Pr-doped ¥28105 (YSO)^. Given the absence 
of atomic motion, extremely long population lifetime can 
be observed. Electromagnetically induced transparency 
and "slow light'— have been demonstrated in rare-earth 
ion-doped crystals, along with storage times greater than 
a second^. 

Among all the rare-earth ions (lanthanides), non- 
Kramers ions with an even number of 4/ electrons exhibit 
the longest coherence lifetime in the absence of external 
magnetic fieldi. With the additional condition on the hy- 
perfine splitting to be of the order of a few tens of MHz 
in the electronic ground state, the rare earth ions are re- 
stricted to praseodymium and europium. Both Pr'^"'" and 
Eu''+ ions have been extensively studied over the past 20 
year o^'^'-'^^i"'^-'^'-'^^'-'^''!"'^^ . However, only dye lasers are avail- 
able at operation wavelengths in Pr'^^ and Eu''+ doped 
crystals. For instance, Pr^^: YSO and Eu"'^: YSO are 
respectively operated at 606 nm and 580 nm. Alterna- 
tively, other non-Kramers rare-earth ions such as thulium 
fall within reach of more tractable lasers. 

Tm^+: Y3AI5O12 (YAG) has been widely studied in 
the field of coherent transient-based signal-processing 
schemeai^ii^ii^iii^ and laserai^ii^. It is especially attrac- 
tive because the ^Hq{Q) — ^Hji{Q) transition at 793 nm 
is in the range of semi-conductor lasers which can be sta- 
bilized to better than 1 kHa^i. 

Although the only isotope of thulium ^^^Tm ions has 
a / = 1/2 nuclear spin, it does not exhibit hyperfine 
structure at zero magnetic field because of J quench- 
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ing22. Lifting the nuclear spin degeneracy with an ex- 
ternal magnetic field offers a way to build a A-system 
provided the magnetic field induces different spin-state 
mixing in ground and excited electronic states. In an ear- 
lier paper— we showed that this was the case for a specific 
magnetic field orientation with respect to the crystalline 
axes. We determined the field direction that gives maxi- 
mum branching ratio between the two optical transition 
probabilities in the A system thus created. In addition, a 
lower bound to the branching ratio was derived from the 
experimental measurement of the gyromagnetic factors. 

Further characterization of this A system includes di- 
rect experimental measurement of the optimal branching 
ratio. This is the purpose of this study. The paper is 
arranged as follows. In Sec. |TT] we summarize the pre- 
vious results concerning the building of a A system in 
Tm"^^: YAG. In Sec. IIIII the experimental setup is de- 
scribed. In Sec. IIVI spectral hole-burning is presented 
and discussed in the case of Tm^+: YAG, and we mea- 
sure the population lifetime of the ground sublevels. In 
Sec. |V]we show that each transition can be individually 
excited and we measure their relative strengths by optical 
nutation and two-pulse photon echoes. From these exper- 
iments we derive an experimental value of the branching 
ratio. 
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II. BUILDING A A SYSTEM IN Tm 



3+. 



YAG 



This section is a brief reminder of the discussion pre- 
sented in a previous paper—. 

Let us consider the optical transition at 793 nm con- 
necting the fundamental levels of the Stark multiplets 
^Hq and ^Hi of Tm^+ in YAG. Application of an ex- 
ternal magnetic field lifts the nuclear spin degeneracy 
by splitting the electronic levels. As shown in Fig. [2l 
a three-level A system would involve the two hyperfine 
sublevels of the ground state ^Hq{Q) and one hyperfine 
sublevcl of the excited state '^H4{Q). If the interaction 
with the magnetic field is restricted to the Zceman effect, 
the electronic levels split into mj = 1/2 and m/ = — 1/2 
spin levels. One of the two optical transitions of the A 
is forbidden. Indeed the selection rule Am/ = forbids 
any electronic excitation to flip the nuclear spin. Fortu- 
nately, the coupling effect of the Zceman electronic and 
the hyperfine interactions mixes the nuclear spin states 
and may lead to comparable transition strengths for the 
two legs of the A if the magnetic field is suitably oriented. 
As a consequence, spin-fiipping transitions |1) — *■ |4) and 
|2) — > |3) are no longer forbidden. The atomic states can 
be represented as the product of an electronic state and 
a nuclear spin state. The ground level nuclear spin states 
are derived from the excited level nuclear spin states by a 
unitary transformation. Therefore the 4 transition dipole 
moments satisfy: 
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FIG. 2: Building a three-level system with nuclear spin levels 
in thulium. 




FIG. 3: Tm substitution sites in the YAG matrix. In each site, 
the transition dipole moment /I is represented by an oblong 
and is directed along the local Oy axis. Here the local frame 
x,y,z is given for site 1. The direction of the magnetic field 
B is defined by angle B. 



The weak (respectively, strong) transition dipole mo- 
ments are referred to as fi^ (respectively, /x^). The 
branching ratio R is the ratio of the transition proba- 
bilities along the two legs of the A: 



R 
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In the YAG crystal, Tm'^"^ ions substitute for yttrium 
ions in 6 crystallographically equivalent but orientation- 
ally inequivalent sites depicted in Fig. O Because of the 
site symmetry (D2), the optical transition dipole moment 
is oriented along the Oy local axis of each site. Let the 
light beam propagate along [110]. Two linear polariza- 
tion directions are considered. 

When it is directed along [111], the optical electric field 
E is perpendicular to the electric dipoles jl of sites 2, 
4 and 6. The only excited sites are then 1, 3 and 5. 
They interact with the light beam with the same Rabi 
frequency f2 = p-'E/h and therefore contribute the same 
amount to the crystal optical density. 
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Alternatively, when the light beam is polarized along 
direction [111], the only excited sites are 1, 4 and 6. Their 
Rabi frequencies are equal. 

For a given magnetic field orientation, the field coordi- 
nates are specific to each site local frame. Since the gyro- 
magnetic tensor is strongly anisotropic in both electronic 
levels, the 6 sites exhibit different hyperfine splittings: 

= v/7^.,eS?+7,%,,i?^+7.Ve5.? (3) 

where g and e denote the ground and excited electronic 
states, and where the gyromagnetic tensor 7 and the 
magnetic field B coordinates are given in the local frame. 
Let B be applied in the bisector plane defined by crys- 
talline axes [001] and [110]. The magnetic field direction 
is defined by angle Q formed with direction [001], as il- 
lustrated in Fig. [3l For any magnetic field orientation 
within the bisector plane, ions in sites 3 and 5 experience 
the same field. Hence they have identical Zeeman split- 
tings and equal branching ratios. Similarly, ions in sites 
4 and 6 experience the same magnetic field and exhibit 
equal branching ratio. 

The branching ratio for ions in sites 3 and 5 was cal- 
culated from theoretical gyromagnetic factors for any 
magnetic field orientation^i. When the magnetic field 
is in the bisector plane. Fig. [4] shows that the theoreti- 
cally predicted branching ratio goes through a maximum 
-Rmax = 0.24. Besides, Fig. U] shows that the accurate 
orientation of the magnetic field is critical around the 
maximum branching ratio orientation. 

The experimental measurement of the gyromagnetic 
tensor coefficients^ made it possible to determine a lower 
boundary for the maximum branching ratio: 

i?max > 0.13 ±0.02 (4) 

These measurements also yielded the magnetic field ori- 
entation that gives maximum branching ratio for ions in 
sites 3 and 5:6 = -49.4 ± 0.2°. This is rather close to 
direction [111] (corresponding to 8 = — arccos(l/-\/3) = 
-54.8°). 

III. EXPERIMENTAL SETUP 
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FIG. 4: Theoretical branching ratio variations for ions in sites 
3 and 5 as a function of magnetic field orientatior>2i. 
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FIG. 5: (color online) Schematic of optical setup after the 
fiber. AO 1, 2 and 3 are acousto-optic devices, S is the 
Tm:YAG sample, and APD is the avalanche photodiode. 

Fig. [5] shows the optical setup after the fiber. Pulse 
sequences are shaped by an Arbitrary Waveform Gener- 
ator (Tektronix AWG 520) that monitors two external 
acousto-optic shifters (AO 1 and AO 2). The acousto- 
optic shifters are imaged on the sample (S) so that the 
beam does not move across the sample when the fre- 
quency is varied. The spot size on the crystal is ad- 
justed to 100 ^m, except in the experiments described 
in Sec. lIVBi where the spot size is 800 fim. The trans- 
mitted light is then collected on an avalanche photodiode 
(APD) HAMAMATSU C5460. An acousto-optic modu- 
lator (AO 3) protects the detector from possible strong 
burning pulses. 



The system is excited and probed with an extended 
cavity semiconductor laser operating at 793 nm. It is 
stabilized on a high-finesse Fabry-Perot cavity through 
a Pound-Drever-Hall servo-loop2i. It is amplified with 
a Toptica BoosTA laser amplifier. A single-mode fiber 
spatially filters the beam. At the end of the fiber one 
measures up to 100 mW output power, depending on the 
amplifier adjustment. 

The 0.1 at.% Tm^+: YAG crystal is 5 mm thick. It is 
cut perpendicular to direction [110] along which the laser 
beam propagates (cf Fig. [3]). We use a commercial Oxford 
6T Spectromag SM4 cryostat. The adjustable magnetic 
field is generated by superconducting coils immersed in 
liquid helium. The sample is cooled down to 1.9 K. 



IV. HOLE BURNING SPECTROSCOPY 

A. Spectral hole burning in an inhomogeneously 
broadened 4-level system 

In thulium, the optical transition at 793 nm exhibits 
a very large inhomogencous broadening Fjnh — 20 GHz. 
Application of a magnetic field splits each electronic level 
into a pair of hyperfine sublevels with a Zeeman splitting 
of the order of 10 to 400 MHz/T, depending on the mag- 
netic field orientation with respect to the local site axes. 
A field lower than 1 Tesla does not modify the absorp- 
tion spectrum since the inhomogencous lincwidth is much 
larger than the Zeeman splitting. Because of inhomogc- 
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FIG. 6: (color online) Optical pumping of four classes of 
ions to the non-resonant ground sublevel by excitation at fre- 
quency fo. Solid lines: strong transitions. Dashed lines: weak 
transitions. Ions are depicted by open circles. 



ncous broadening, a laser excitation at given frequency 
i/Q simultaneously excites 4 different classes of ions, as il- 
lustrated in Fig. [6l Each class refers to a specific ground 
and excited state sublevels that are connected by the in- 
coming field. Initially, in each ion class both ground sub- 
levels are equally populated. This is even true at working 
temperature (T = 2 K), since the level splitting is much 
smaller than ksT/h ~ 30 GHz. 

When the sample is exposed to sustained incident light 
at I'o, it experiences optical pumping that tends to unbal- 
ance the ground level state distribution. Ions in the reso- 
nantly excited sublevel are transferred to the off-resonant 
sublevel (cf Fig. [6]). This gives rise to a hole in the ab- 
sorption spectrum at frequency uq. Besides, since we 
consider 4-level systems, additional features appear. Two 
side holes are observed at I'o ± Ae, where the probe field 
excites transitions from the depleted ground sublevel to 
the other excited state sublevel. Conversely, transitions 
from an overpopulated ground sublevel to either excited 
sublevel give rise to antiholes in the absorption spectrum. 
Such antiholes are observed at VQ±/S.g and i'o±(Ac,± Ag). 

Depending on the sublevels they connect, the various 
transitions excited by the probe beam can be sorted as 
weak or strong ones. It should be noted that at least 
one strong transition contributes to any hole and anti- 
hole, except for the antiholes located at vq ± (Ag -I- Ae). 
These antiholes arise from a probing process along a weak 
transition in the ion classes (ii) and (iv) of Fig. [51 As a 
result, the observation of an antihole at this position shall 
bear evidence that the nuclear spin flip selection rule is 
actually relaxed. 

We perform spectral hole-burning in Tnr^+iYAG. A 
0.19 T magnetic field is oriented in the direction that 
is expected to give maximum branching ratio for ions 
in sites 3 and 5 (6 = -49.4°, cf Sec. |Tl|. The laser 
beam is linearly polarized in direction [111] so that only 
sites 1, 3 and 5 are excited. The burning step consists 
in a series of ten 50 ^s pulses every 10 ms. Then the 
sample is probed with a weak 1 ms long pulse chirped over 
a 35 MHz interval. The burning and probing sequence 
is repeated every 150 ms. Transmitted probe intensity 
is collected on a photodiodc, averaged over 8 shots and 
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FIG. 7: Hole-burning transmission spectrum for Tm"^"*": YAG 
under a 0.19 T magnetic field. The spectral holes and an- 
tiholes correspond to site 3 and 5 ions. The features corre- 
sponding to site 1 are out of the probing spectral window. 



plotted in Fig. [T] The holes and antiholes are ascribed 
to sites 3 and 5. We measured the ground and excited 
splittings for these two sites in this field direction: 



Ag/B = 38.2 ± 1 MHz/T 
Ae/S = 15.5 ±0.7 MHz/T 



(5) 
(6) 



where B is the magnetic field amplitude. This agrees 
with the calculated splittings derived from the exper- 
imental gyromagnetic tensor components'^: Ag/S = 
36.0 ± 1 MHz/T and Ae/B = 16.0 ± 0.8 MHz/T. The 
ground level gyromagnetic tensor values are close to those 
measured by Schmidt in TmAG by Nuclear Magnetic 
Resonance2^. 

One can see in Fig. [7] that the antiholes arc generally 
wider than the holes. This suggests that the inhomoge- 
neous broadening of the Zeeman splitting is larger in the 
ground state than in the excited state. 

The hole-burning features associated with site 1 cannot 
be seen in Fig. [7] because the Zeeman splitting for site 1 
is over 60 MHz/T for this magnetic field orientation. 



B. Sublevel population lifetime 

The hole-burning spectrum is erased by return to ther- 
mal equilibrium. To measure the corresponding decay 
rate, we measure the central hole depth as a function of 
delay between pumping and probing steps. The sample 
is cooled down to T = 1.9 K. A 0.45 T magnetic field 
is applied in direction 8 = —54.5°, very close to [111]. 
It should be noted that this is not the optimal direction 
discussed in Sec. [Hi 

The light beam is polarized in direction [111] so that 
only sites 1, 4 and 6 are excited, with identical Rabi fre- 



5 



e 



-54.5° e = -49.4° (maximum R) 



1 

0,8 
0,6 



0,2 



2 4 6 8 10 

Time delay between burning and probing (s) 



FIG. 8: Decay of central spectral hole after burning operation 
for sites 1, 4 and 6 (triangles) and sites 3 and 5 (circles). 
For these two measurements an exponential function is fitted 
(dashed lines). In this experiment only, the 0.45 T magnetic 
field is along [111]. 



site 1 328 MHz/T 
sites 3, 5 15 MHz/T 



306 MHz/T 
36 MHz/T 



TABLE I: Calculated values for normalized ground splitting 
Ag/B in different sites, for two specific field orientations. 



e = -54.5° to e = -49.4°, the branching ratio of the 
optical transitions strongly increases (cf Fig. 2]) , but the 
level splitting Ag varies little, as shown in Table [J The 
decay rate is expected to be mainly sensitive to this level 
splitting. 

Before sites 3 and 5 can be coherently manipulated for 
quantum storage, it is necessary to prepare site 1 ions 
so that they do not interfere. Their very long popu- 
lation lifetime guarantees that a bleaching process can 
store them for a long time in a non-resonant ground sub- 
level. As for sites 3 and 5, a spin level lifetime of a few 
seconds is long enough for coherent manipulation. 



V. BRANCHING RATIO MEASUREMENT 



quency. In addition, since they are magnetically equiva- 
lent, their Zccman splittings are equal (A^ = 181.4 MHz 
and Ae = 36.9 MHz) and they exhibit the same popula- 
tion lifetime. Fig. [8] shows that the central hole docs not 
decay for 10 seconds following excitation. In other words, 
the population lifetime relative to sites 1, 4 and 6 for this 
magnetic field orientation and strength is estimated to a 
few minutes. 

The light beam is then polarized along [111] so as to ex- 
cite sites 1, 3 and 5 with identical Rabi frequency. Since 
we are interested in measuring the sublevel population 
lifetime for sites 3 and 5 only, we isolate ions in sites 3 
and 5 by means of a bleaching process: a strong beam 
chirped over a 25 MHz interval optically pumps a large 
spectrum range of ions in site 1 to the non-resonant very 
long lifetime ground sublevel. Site 3 and 5 ions also un- 
dergo optical pumping but the chirp amplitude is much 
larger than the ground sublevel splitting Ag = 6.9 MHz. 
As a consequence, a large fraction of ions in sites 3 and 5 
is pumped back to the initial sublevel. Three seconds af- 
ter this bleaching process, we perform ordinary spectral 
hole-burning and observe the decay of the central hole 
depth by varying the delay between burning and probing 
(cf Fig. [8]). An exponential fit gives a population lifetime 
of 4.5 ± 0.5 s for sites 3 and 5. 

Therefore in crystallographically equivalent sites 1,3 
and 5, the sublevel population lifetime apparently in- 
creases with the ground state splitting. At present we 
cannot explain this behavior. In contradiction with our 
measurements, N. Ohlsson et al. observed longer life- 
times for smaller splittings in a given site^^. 

Since the spin level thermalization does not involve op- 
tical transitions, we expect the corresponding decay rate 
to be insensitive to a small angle rotation of the mag- 
netic field. Indeed, when B orientation is varied from 



In order to precisely compare the optical transitions 
along the two legs of the A, we resort to coherent pro- 
cesses along with frequency selective optical pumping. 
This enables us to isolate the weak optical transition con- 
tribution. 



A. Optical nutation 

A two-level system driven by a resonant monochro- 
matic field E oscillates between its ground and excited 
states at Rabi frequency fl ~ jl ■ E/h. Along with this 
population oscillation, the system emits an induced ra- 
diation that modulates the driving field at frequency Q. 
This gives a direct and robust way to measure the Rabi 
frequency of a transition even when the inhomogeneous 
broadening is infinite and the beam spatial profile is not 
uniform: for an infinitely inhomogeneously broadened 
medium driven by a stepfunction gaussian laser beam at 
t = 0, the transmitted signal takes the form of damped 
oscillations21: 

:^ = 1-2(1-10-^):MM fori>0 (7) 

where Iq is the incoming beam intensity, D is the opti- 
cal density of the sample, il is the Rabi frequency of the 
ion-laser interaction in the center of the beam and Ji is 
the Bessel function of order 1. This expression is valid 
if the optical density of the sample is small {D < 0.5). 
According to Eq. [71 the first maximum of /nut (t) is lo- 
cated at tmax = 5.1/f2, therefore the Rabi frequency can 
be directly obtained from the nutation signal. 

The optical density D can be assessed by comparing 
the transmission of the material at i = with its trans- 
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FIG. 9: Optical nutation signal with zero magnetic field (open 
circles) when a step-function monochromatic excitation is ap- 
plied at t = 0. The signal is fitted according to Eq. [7] with 
Q = 2.6 MHz and D = 0.32 (solid line). 



mission after a long exposure to laser light: 

/nut(i OO) 



D = log 



10 



/nut(t-O) 



(8) 



This way, one can measure the optical density of the ma- 
terial with a single nutation temporal profile. Besides, 
the optical density is proportional to and to the num- 
ber of active atoms. If optical nutation is performed on 
each transition, then the optical density measurement 
gives direct access to the relative dipole moments of the 
weak and strong transitions. 

We first observe optical nutation on the two-level sys- 
tem formed by Tm"^^: YAG at zero magnetic field. The 
transmitted signal is displayed in Fig. [9l The sample is 
cooled down to T = 3 K and is excited with a 10 /iS res- 
onant light pulse polarized along [111]. For this light po- 
larization, sites 1 , 3 and 5 participate in absorption with 
equal Rabi frequency and each one of them contributes 
one third to the total optical density Dq = 0.32. 

Let us apply a magnetic field in the direction that is 
expected to give maximum branching ratio for sites 3 
and 5 (cf Sec. |TI]). Application of a magnetic field splits 
the electronic levels, the ions being equally distributed in 
both ground sublevels. As mentioned in Sec. IIV Al this 
does not change the absorption profile of the sample, so 
the nutation signal is not modified. 

We now apply 10 /is pulses. The sample temperature is 
decreased gradually from 5 K to 2 K. As long as thermal- 
ization processes dominate, the system returns to thermal 
equilibrium between successive pulses. As the tempera- 
ture goes down, thermalization gets slower. Ultimately, 
resonant sublevel depletion by optical pumping super- 
sedes repopulation by relaxation processes (see Fig. [6]) 
and the nutation signal disappears. 

The repopulation process can be assisted by opti- 
cal repumping. If one tunes the repumping laser to 



i^R = t^o ~ so as to counterbalance optical pump- 
ing at vqi only ions in sites 3 and 5 from classes (i) and 
(ii) are pumped back into the initial ground sublevel [see 
Fig. [TUTa)]. All ions in site 1, as well as ions in sites 3 
and 5 with excitation schemes (iii) and (iv) , stay off reso- 
nance with the repumping beam. In a given class of ions, 
the steady-state fraction of ions in the resonant sublevel 
reads as 



(9) 



P 



2k 



where p, r and k respectively stand for the pumping rate, 
the repumping rate and the temperature dependent re- 
laxation rate. In the absence of optical excitation this 
fraction reduces to p = 1/2, which means that the ground 
state population is evenly shared between the two sub- 
levels. When r = and p ^ k, all ions are pumped to 
the off- resonance sublevel and p = 0. When r,p^ k the 
relative population reduces to p = r/(p + r). 

Optical pumping consists in repeated excitation and 
relaxation processes. Excitation is carried out by a se- 
ries of 10 fis pulses. Within this 10 p,s pulse duration the 
transition at vq is bleached, resulting in evenly populated 
ground and upper levels. As will be shown later in this 
section, this is true even when the excitation takes place 
along the weak transition. Therefore the excitation pro- 
cess brings the system in the same final state, whatever 
excitation path is followed. 

Relaxation to the ground state strongly differs from 
spontaneous emission decay. Indeed, non-radiative pro- 
cesses dominate, proceeding through the intermediate 
states ^H^ and ^F4. We assume here that these step- 
wise relaxation processes relax the nuclear spin selection 
rule. To sum up, neither excitation nor relaxation de- 
pends on the transition strengths, and the same value of 
p is used for classes (i) and (ii). 

The electronic transition dipole moment docs not de- 
part from its zero magnetic field value. Since only 2 sites 
out of 3 participate in the absorption, the optical density 
at stake here reads as 



Ds+w — 



2 



(10) 

where Dq is the optical density at zero magnetic field. 

Similarly, if the repumping laser is tuned to h'R = h'a + 
Ag + Ae [cf Fig. [TUTb)]. only ions from system (iv) are 
repumped to the resonant sublevel: the nutation signal 
is ascribed to ions in sites 3 and 5 resonating on a weak 
transition. The optical density is expected to be much 
smaller than Dq and reads as 



Du, = -:Dq 



2 

3^°r 



R 



(11) 



Finally, if i^ji = vq + Ag — Ag, only ions from system 
(iii) are repumped to the resonant sublevel: the nutation 
signal is ascribed to ions in sites 3 and 5 resonating on a 
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FIG. 10: (color online) Examples of repumping schemes. Four 
classes of ions are excited at frequency vq. The repumping 
frequencies are displayed as thick solid (respectively, dashed) 
arrows along strong (respectively, weak) transitions, (a) ur — 
vq — Ag-. only systems (i) and (ii) participate in the nutation 
signal, along a strong transition for system (i) and along a 
weak transition for system (ii). (b) vr — vo + Ag + A^: only 
ions in system (iv) participate in the nutation signal, along a 
weak transition. 



strong transition. The optical density reads as 



2 1 
-Do 

3 °1 + 



R 



(12) 
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FIG. 11: (color online) Optical nutation signals for various 
repumping frequencies. The four solid traces are obtained 
with the same incoming laser intensity, a: with zero magnetic 
field, b: with magnetic field and repumping at frequency i^r — 
vq — Ag, c: with magnetic field and repumping at frequency 
vr = I'o — Ag — Ae, d: with magnetic field but no repumping 
beam. The dashed traces are theoretical curves given by Eq.[7] 
and fitted to the experimental data. The graphs are vertically 
offset for clarity. 



The combination of optical pumping and repumping 
processes makes it possible to isolate a nutation signal 
related either to a weak or to a strong transition. The 
transmitted intensity for various repumping frequency 
values is displayed in Fig. [TlJ The repetition rate is 
50 s~^, the sample temperature is T = 2 K and the 
incident intensity is the same whatever the repumping 
frequency value. 

The optical density of the sample under zero mag- 
netic field is Dq = 0.36 according to Eq. [S] applied to 
Fig. [TlTa). When a 0.3 T magnetic field is applied, the 
sample becomes transparent: the damped oscillations 
disappear and the transmitted intensity is almost con- 
stant [cf Fig. md)] . 

The repumping process consists in a series of ten 100 /is 
chirps over 1 MHz around the repumping frequency ve.- 
It is applied between successive excitation pulses. In 
Fig. fTlTb) we plotted the transmitted intensity when 
a repumping beam is applied at frequency vq — Ag. 
The optical density is measured on the nutation signal: 
Ds+w = 0.167. Eq. [TO] gives the relative population of 
the resonant sublevel in systems (i) and (ii): p = 69%. 
Besides, in both cases (a) and (b) of Fig. [Til the nuta- 
tion signal is associated with exciting strong transitions 
as well as weak transitions. This is why the Rabi fre- 
quency and hence the position of the first maximum in 
the nutation signal do not change. 

Finally, in Fig.fTlTc) we plotted the transmitted inten- 
sity for a repumping beam at i/q + Ag + Ag. The optical 
density of the sample excited only along a weak transi- 
tion is very small: D^j = 0.018 ± 0.004. No damped os- 
cillations are visible, but the transmitted intensity slowly 
increases with time during excitation. The experimental 



data clearly show that the position of the first maximum 
in nutation signal (c) is different from cases (a) and (b), 
even though it cannot be accurately located. This proves 
that signal (c) does not come from a residual signal emit- 
ted along a strong transition. 

According to Egs.fTOlandfTTl the ratio of Z?^ and Ds+w 
finally leads to: 

u 2 

i? = 0.12 ±0.03. (13) 

Ms 

which is consistent with the lower boundary given in Eq.[4| 
derived from gyromagnetic tensor measurements^^. 

The position of the first maximum of nutation profile 
for the weak transition is -v/i? times the corresponding 
position for the strong transition. According to the mea- 
sured R value, this maximum occurs at i ~ 5.5 ^s. There- 
fore, as previously assumed, saturation is reached within 
the 10 fis pulse duration even for the weak transition. 

B. Photon echo 

In this section, two-pulse photon echoes are achieved 
separately on each leg of the A. Application of a magnetic 
field splits the electronic levels, and repeated excitation 
achieves optical pumping between the two ground sub- 
levels. Therefore the echo is present as long as the tem- 
perature is high enough so that relaxation compensates 
for optical pumping. 

Fig. [T2] shows a two-pulse photon echo sequence in 
Tm:YAG under 0.3 T magnetic field, with a 50 s~^ rep- 
etition rate. The two pulses are gaussian and have equal 
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FIG. 12; Photon echo in Tm='+: YAG under magnetic field 
with two gaussian pulses at T = 2.6 K (solid line) and T — 
1.9 K (dashed line). The detection channel is gated by an 
acousto-optic modulator (AO 3 in Fig. [5]) which attenuates the 
excitation pulses (grey area). This AOM opens a few fis after 
the second pulse to let the photon echo through. The echo 
disappears when optical pumping dominates thermalization. 



intensity. At T = 2.6 K, thermalization still dominates 
the optical pumping caused by repeated excitation. In 
each two-pulse sequence, the first pulse partly bleaches 
the optical transition, so the second pulse is less absorbed 
than the first one. As the sample is cooled down, optical 
pumping becomes dominant. At T = 1.9 K, the reso- 
nant level is totally empty and the photon echo signal 
disappears. The excitation pulses propagate through a 
transparent medium where there is no bleaching process, 
which is why their transmitted intensities are equal. 

We need to perform the same repumping process as 
described in Sec. IV Al in order to recover a signal and 
ascribe it to a weak or strong transition. When the re- 
pumping frequency is tuned to vq ± (Ag — Ag), ions in 
system (i) or (iii) in Fig. [TO] are repumped and the pho- 
ton echo is emitted along the strong transition at fre- 
quency Vq. When the repumping frequency is tuned to 
i/Q ± (Ag -I- Ae), the photon echo is emitted along the 
weak transition at frequency vq in ions in system (ii) or 
(iv). To make sure that the detected signal is imputable 
to the weak transition, we slightly detune the repumping 
beam: the echo intensity goes down, which proves that 
the off-resonant repumping is less efficient and that the 
signal does not come from a residual signal emitted along 
a strong transition. 

The photon echo amplitude can be expressed as a prod- 
uct of the transition dipole moment and a function of 
the Rabi frequency. This applies to both transitions, ei- 
ther strong or weak. Since the Rabi frequency is propor- 
tional to the transition dipole moment and to the applied 
field, the photon echo intensity along the strong transi- 
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FIG. 13: Weighted echo intensity on the strong transition for 
difi'erent excitation pulse intensities (squares). Theses values 
are interpolated with a straight line (solid line) in order to 
solve Eq. 1171 A graphical solution is given by the intersection 
point of W(T) = J^X) xJ/Jo and W{I) = I^{To). The value 
is directly read out as -R = I/Tq at the intersection. 

tion reads as: 

^ ^il g{^ilI) (14) 

where T represents the laser intensity and where the 
shape of function g need not be detailed any further. 

Conversely, when the weak transition is excited with 
the same intensity X, the photon echo intensity reads as 

i^{T) = fil gifili) = i? gip^'sR^) = R UR^) (15) 

where R is the branching ratio. 

The photon echo intensity 2s on the strong transition is 
measured for different excitation intensities. A weighted 
echo intensity is defined as: 

M^(X) =X,(J) xX/Xo (16) 

This function is interpolated with a polynomial function 
as shown in Fig. 1131 If we measure X^, on the weak tran- 
sition for a given intensity Xq, wc arc able to derive the 
branching ratio value by solving the following equation: 

X^(Xo) = i?X,(i?Xo) (17) 

As shown in Fig. [13] and according to Eqs. [16] and [TTl 
the graphical solution of this equation is given by the 
intersection point of W{X) = Xs(X) x X/Xq and W{1) ~ 
X„(Xo). The x-coordinatc of this intersection point is the 
branching ratio R: 

i? = 0.130 ±0.015. (18) 

This is consistent too with the expectations derived from 
gyromagnetic factor measurements^^. 



Discussion VI. CONCLUSION 



With optical nutation and photon echo experiments we 
measured two branching ratio values that are in agree- 
ment with one another. 

As far as theoretical gyromagnetic tensors are con- 
cerned, the {x, y) anisotropy is much larger in the ground 
state than in the excited stated. This difference, known 
as anisotropy disparity, is at the origin of the high value 
for the maximum theoretical branching ratio -Rmax = 
0.24. The experimental value of anisotropy disparity ap- 
peared to be smaller than predicted by theory^^. As a 
result, a branching ratio value close to its lower boundary 
0.13 was expected. The direct measurements presented 
in this paper confirm this result. 

These measurements have been achieved for the mag- 
netic field orientation (6 = 49.4 ±0.2°) given by de Seze 
et al^. We did not study the dependence of R as a 
function of O. However, the branching ratio is almost 
constant over a 2° wide interval^i. Therefore, even with 
a 0.5° mechanical precision on the angle, one can reason- 
ably consider that the measured branching ratio is close 
to its maximum. 

Finally, in the experiments presented in Sec. IV Ai the 
nutation signal along the strong transition goes through 
its first maximum at t = 2 fis. This means that a 7r-pulse 
along the strong transition would take 1.2 fis. According 
to the present branching ratio measurements, a 7r-pulse 
along the weak transition would then take 1.2/^/R = 
3.4 ^s. Therefore, we are able to excite each leg of the 
A system with a 7r-pulse, which is a capital asset in the 
perspective of coherent processes dynamics. 



We directly measure the transition probability ratio 
along the two legs of the A through coherent processes 
along with frequency selective optical pumping. This is 
performed in the optimum magnetic field orientation de- 
rived from a previous experimeni*^^. The direct mea- 
surement of the branching ratio is consistent with the 
value derived from the gyromagnetic tensor components 
and confirms that a A-system can operate efficiently in 
Tm:YAG. This way we have cross-checked and supported 
the theoretical description, complementing prior experi- 
mental data^ with a new set of data different in nature. 
In addition, this provides us with the branching ratio 
value we have to know in order to efficiently build and 
control nuclear magnetic superposition states with opti- 
cal excitation. 

Side investigation on the sublevel lifetime leads to un- 
expected results. Indeed, broadly spaced sublevels ap- 
pear to be more stable than closely spaced ones, in con- 
tradiction with previously reported experiments^. This 
has to be clarified. Still, the measured lifetime is consis- 
tent with long-time coherent manipulation of spin states. 
The demonstrated existence of the A system promises in- 
teresting spin state properties in either ground or excited 
level. Next step will deal with the investigation of these 
coherent spin features. 
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